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Abstract: In this paper, a liquid multi-parameter decoupling method with only one Lamb 
wave sensor is presented. In a Lamb wave sensor, antisymmetric modes (A 0 i mode for low 
frequency, A 0 3 mode for high frequency) and symmetric modes (So mode) are used to 
detect multiple parameters of a liquid, such as its density, sound velocity, and viscosity. 
We found they can play very different roles in the detections. For example, the Aoi mode is 
very sensitive to the liquid's density but the A03 mode is sensitive to the sound velocity. 
Here, the A 0 mode is used to identify the density of the detected liquid and with this 
density value we obtained the viscosity by the amplitude shifts of the So mode. This could 
be a way to distinguish an unknown liquid with high sensitivity or to solve the problem of 
selectivity of label-free detection on biosensors. 

Keywords: Lamb wave sensor; density; sound velocity; viscosity; determination of 
multi-parameters 
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1. Introduction 

Acoustic sensors have been widely used in chemical/biological fields with the label-free detection 
method to detect the mass changes on the sensor surface [1-17]. However, the selectivity of this 
method is often poor due to the absorption of non-target molecules, which are difficult to distinguish 
using only one -parameter sensors [18,19]. This has been a common problem of all kinds of label-free 
biosensors. The detection of multiple parameters using a multi-mode acoustic sensor significantly 
improves the label-free detection method. Acoustic waves which travel in a medium can have multiple 
modes and this character has been successfully used in a combined detection of density and viscosity 
for high viscosity solutions [20]. However, for low viscosity solutions like aqueous electrolytic 
solutions or bio liquids, it is still a challenge [9,20-24]. 

The micro Lamb wave sensor is a powerful tool for liquid detection because it is easier to get 
multi-mode vibration, and it performs with high sensitivity and low attenuation [25-29]. The two 
well-known basic modes, the antisymmetric mode (A 0 mode) and the symmetric mode (So mode), have 
already been used for solving the problem of temperature compensation on a chip by the authors 
of [26,30]. The A 0 mode has been used to measure the concentration of bio/chemical liquids, such as 
the concentration of methanol for direct methanol cell applications [18]. However, it shows uncertain 
frequency shifts direction with the concentration changes of the liquid [18,31]. In fact, the parameters 
of a liquid, like the density, acoustic sound velocity and viscosity, will work together for a mode. So 
it's hard to decouple them with only one mode. 

In this work, both the antisymmetric mode and the symmetric modes are used to decouple the liquid 
physical parameters. The characters of A03 mode, with a wavelength of about one-third of A 0 i mode 
(low frequency of A 0 mode), are investigated with the liquid loading for the first time. In order to 
discuss the response of Lamb wave sensor to the liquid loading, two types of experiments were set up. 
The first experiment was used to measure the frequency shifts with the same type of solutions with 
different concentrations loaded on the sensors' surface, such as the solutions of NaCl. It shows that the 
characters of the modes are very different. Although both A01 and A03 modes belong to Ao mode, their 
frequency shifts showed an opposite behavior with the increase of the concentration of the loading 
solutions. The second experiment was used to measure the frequency changes with different types of 
solutions and different concentrations, such as solutions of KC1, NaBr and KBr. Although the tested 
objects are two different kinds of solutions with different concentrations, the frequency of A01 mode 
shows almost has no movement and A03 mode frequency shows a great difference. The frequency 
shifts of So mode for these solutions are not apparent. These will make it possible to be only sensitive 
to the liquid viscosity except for the environmental temperature. This essay attempts to decouple the 
functions of A 0 i, A 0 3 and S 0 modes to obtain the basic physical parameters (density and sound velocity, 
viscosity) with different liquid loading responses. From the values which have been reported from 
literature [32] we can also determine the types of the solution. This work provides a selective, sensitive 
method for the measurement of physical parameters to an unknown solution. It can be useful for 
label-free biosensors. 
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2. Working Principle 

Multi-mode wave, such as A 0 i mode, So mode and their harmonic ones (A 0 3 mode, for example) can 
be excited in Lamb wave sensors and detected directly by using a pair of inter-digital transducers 
(IDTs) [33] located on the surface of the piezoelectric layer [34,35]. In liquid sensing, the liquid and 
IDT will be located on the two opposite sides of the membrane. For the Aoi mode and the A03 mode, 
when the Lamb wave phase velocity is less than the liquid sound velocity, there are evanescent waves 
produced around the membrane-liquid interface. The phase velocity of the Lamb waves within the 
evanescent penetration field will be the same with the phase velocity inside the membrane. Taking 
account of the bending stiffness and the in-plane tension of the plate (B;) for the A 0 i mode and the A03 
mode, the phase velocity will be [31,36]: 

c Pi = pJ(M + m Ll ) (1) 

where i denotes 1 or 3 for the A 0 ; mode, B\ reflects the influences of the bending stiffness and the 
in-plane tension of the plate for the Ao; mode, M is the mass per unit area of the plate, the effective 
mass m Li in the so-called evanescent penetration depth S E i equals: 

™u = Pl^Ei (2) 

In which, p L is the detected liquid density. The evanescent penetration depth S E i obeys the following 
equation [37-39]: 

S EI =l/(2xJl-(c„/c L ) 2 ) (3) 

where A; and c Pi denote the Lamb wave wavelength and phase velocity respectively in each mode. c L is 
the bulk acoustic velocity of the detected liquid. c Pi is determined by the frequency (fi) via c Pi =f t A,-. 
By substituting the Equations (2) and (3) into the Equation (1), we will get the following formula: 

pA/(2^l-(c n /c L ) 2 ) = Bjc 2 Pl -M (4) 

In this formula, the constants B h A,, and M are independent of the liquid type; A, is determined by the 
structure of the device. With the implicit Equation (4), the density and sound velocity of the liquid on 
the sensor can be decoupled when the frequency response of the A 0 i mode and the A03 mode are 
obtained by the experiments. 

In the measurement, the relative frequency shifts (A/7/) will be used frequently. By taking the term 
cpi =fi A; into Equation (1) and deriving of the equation, the value of A/7/) for the Aoi mode and the A03 
mode will be: 

A/7/; = -A(p L S Ei )/2M = -(SeAPl + p L M a )/2M (5) 

Obviously, the density and the sound velocity are the two factors affecting the relative frequency 
shifts. When the phase velocity is far less than the liquid sound velocity (cp « c£), such as the Aoi 
mode, the values of Se and ASe approach A/27T and 0, respectively. Then, in this situation the variations 
of the liquid sound velocity will have an almost negligible influence on the relative frequency shift. 
When the phase velocity is close to the liquid sound velocity, such as the A03 mode, the value 
1 - (cp/c L ) 2 approaches 0. The influence of the terms S E and c L cannot be neglected any more. The 
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liquid density and sound velocity will simultaneously affect the frequency shifts. The value of Aflf can 
be positive or negative which depends on the sum of the value of plASe and the value of Se^Pl- 
Anyway, if we combine these two cases (A 0 i mode, A 0 3 mode) simultaneously, we can get the values 
p L and c L with high sensitivity. 

For the So mode, when the wavelength of the Lamb mode is larger than the thickness of the plate, 
the phase velocity (c p ) for the principal symmetric mode can be simplified as [40]: 



c p =c, 
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where q and c t are the velocities of longitudinal (dilatational) and transverse (shear) waves of the 
membrane, c L is the liquid sound velocity, p m and p L are the density of the solid membrane and liquid 
respectively, d is the membrane thickness, / denotes the square root of — 1. 

Evidently, the effect of liquid on the propagation of the So mode does not change the real part of the 
phase velocity, but adds a very small attenuation of the amplitude [40]. This mode is suitable for the 
detection of attenuation. The amplitude (A L ) response of an unknown solution can be expressed by: 

A L =\e~ a ^ (7) 



where the attenuation coefficient o. X l is proportional to (pi?//,) 0 5 , with n L being the viscosity of the 
liquid. Similarly, the amplitude response (A w ) of water is given by: 

A w =A 0 e~ a ^ (8) 

In engineering, the insertion loss (dB) is widely used. The values of A L and A w can be transformed 
into values in dB scale which are denoted by Aub and A W dB, respectively. Therefore, the attenuation 
difference (AAjb) between an unknown solution and water can be expressed in dB scale as follows: 

Ms = A idB ~ A wdB = Yx {PlVl f 2 + 72 (9) 

in which both of the slope y\ and yi are constant, and the constant 72 is determined by the reference 
liquid (water). As the density and the sound velocity of the solution can be estimated by both the 
frequency shifts of the A 0 i mode and the A 0 3 mode, the multi-parameters (pi, cl, Ml) of the solution 
can be decoupled simultaneously with these three modes. According to these principles, a series 
of experiments were set up to investigate the responses of the multi modes of Lamb wave to the 
loading solutions. 

3. Experiments for Liquid Detection 

3.1. Experimental Setup 

The micro Lamb wave device in Figure 1(a) contains a silicon membrane (length 7.8 mm, thickness 
-12 (am) with a ground layer (Ti/Mo, GND, -0.2 jam) and a piezoelectric layer (aluminum nitride, 
A1N, -1.8 (am). Lamb waves are excited and detected directly using bidirectional inter-digital 
transducers (IDTs) [33] located on the surface of the A1N layer. There are six pairs of fingers on 
bidirectional IDTs in each exciting and detecting transducer. The period of bidirectional IDT is about 
400 jam. As waves are partly reflected at the end of length-limited membrane (-7.8 mm), the device 
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has strong signal without reflectors on both ends of the membrane. With the layers described as 
in [34], the mass per unit area of the membrane (M) is about 0.0355 kg/m 2 . 

The micro Lamb wave device is packaged directly with the printed circuit board (PCB), as shown in 
Figure 1(b). The network analyzer (Agilent 4395 A), connected with the PCB, is used to excite and 
receive the acoustic signals. The device is protected with one cover on top of the system. The tube is 
used to pass the liquid into/out of the chamber, which is sealed up with the polymethyl methacrylate 
(PMMA) cover under the PCB. 

Figure 1. The system for liquid detection, (a) Schematic diagram of the micro Lamb wave 
sensor interaction with liquid, pc density, cl. sound velocity, rji- viscosity; (b) Micro Lamb 
wave sensor packaged with printed circuit board (PCB). 




Figure 2. The responses of the Aoi mode, the A 0 3 mode, and the So mode to the air and the 
water loading on the Lamb wave sensor. 




Frequency (MHz) 



When the device is loaded with air or water, multi-modes can be excited and detected effectively, 
including the Aoi mode, the A 0 3 mode, and the So mode (Figure 2). As the Aoi and A 0 3 modes are 
harmonic, the wavelength of the A03 mode is about one third of the wavelength of the Aoi mode. 
According to the positions of the A 0 i mode, A 0 3 mode and S 0 mode in the dispersion curves [41,42], the 
resonant frequencies if) of these three modes should have the following relation: f A oi <fA03 < fso- The 
measured resonant frequency of the A 0 i and A03 modes for water loading are 0.987 and 9.233 MHz, 
respectively. With the frequency if) and the wavelength (X), the corresponding phase velocities (fit) of 
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these two modes are 385 m/s and 1,200 m/s respectively. Although the phase velocity of the Aoi mode is 
far less than the sound velocity of water (1,485.5 m/s), it is close to the sound velocity of water for the 
Ao3 mode. With the same water loading for the So mode, the measured central resonant frequency is 
about 20.86 MHz. The corresponding phase velocity of the S 0 mode is about 8,135.4 m/s. Further 
experiments were conducted to investigate the response of multi-modes to the different solutions. 

3.2. Experiments for One Species Solutions (NaCl Solutions) with Different Concentrations 

The first experiment was done to measure known solutions with different concentrations (Figure 3), 
such as NaCl solutions. When water is designated as the reference liquid, the relative frequency shifts 
(Af/f = (Solution - /water)//water) with concentration are different for these three modes (Aoi mode, 
A03 mode, and S 0 mode), where / so i u tion and f,^ are the measured frequencies for the solution and the 
water, respectively. The frequency of the Aoi mode decreases with the concentrations and the A 0 3 
mode behaves oppositely. In the case of the A 0 3 mode of NaCl solutions measured in Figure 3, the 
absolute value pASg is bigger than the absolute value SgAp (Table 1). This causes the positive 
frequency shifts of A03 mode in NaCl solutions measurements. In any case, when the frequency shift is 
positive, this means the Lamb wave phase velocity is close to the liquid sound velocity. In this case, 
this phenomenon can be used to decouple the liquid sound velocity. 

Figure 3. Relative frequency shifts (Af/f) for the Aoi mode, the A03 mode, and the So mode 
in the measurements of the NaCl solutions. 




0 5 10 15 20 

NaCl mass percent in solutions (%) 



Table 1. Based on the measured frequency shifts of NaCl solutions measured in Figure 3, 
comparison of the value pAdg and the value dgAp for the A 0 3 mode. 



Concentration (%) 


Density (10 3 kg/m 3 ) 


c P (m/s) 


S E (fim) 


S E Ap 


pAS E 


0 


0.9981 


1,200.30 


35.12 


0 


0 


0.99 


1.0052 


1,201.71 


34.70 


0.25 


-0.42 


5.67 


1.0389 


1,206.55 


32.81 


1.34 


-2.39 


9.09 


1.0640 


1,209.03 


31.58 


2.08 


-3.77 


13.79 


1.0995 


1,211.59 


30.01 


3.04 


-5.62 


16.64 


1.1212 


1,212.69 


29.13 


3.59 


-6.71 


19.33 


1.1425 


1,213.53 


28.40 


4.10 


-7.67 



Compared with the frequency shifts measured with the Aoi mode and A03 mode, the frequency of 
the So mode shows negligible shift with different concentrations (Figure 3). 
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3.3. Experiments for Three Different Unknown Species Solutions 

The second experiment was done to measure the modes' responses to three different species of 
aqueous electrolytic solutions with different concentrations (Table 2). We will check the possibility of 
the method to decouple the density and the acoustic sound velocity of these solutions with the 
measurements of A 0 i mode and A03 mode. Three different species of aqueous electrolytic solutions are 
prepared, they are listed as solution A, solution B, and solution C respectively, as it is shown in Table 2. 
From No. 1 to No. 6 solutions, they are the solution A. No. 7 and No. 8 solutions are the solution B. 
No. 9 and No. 10 are solution C. In each kind of solution, the concentration increases with the serial 
number. The So mode is still insensitive to the solution changes compared with the measurement in 
water, and its values are not listed here. 



Table 2. Relative frequency shifts (Af/f) for the A 0 i mode and the A03 mode in the 
measurements of three kinds of unknown solutions. In each kind solution, the concentration 
increases with the serial number. 









Solution A 






Solution B 


Solution C 


No. of solutions 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Concentration (%) 


0.99 


6.05 


9.97 


16.03 


19.96 


23.95 


9.35 


18.07 


8.98 


17.32 


Aflf (A01 mode) 
(-1 x 10 4 ppm) 


0.236 


1.09 


1.846 


2.886 


3.654 


4.434 


2.243 


4.474 


2.216 


4.356 


Aflf (A03 mode) 
(10 4 ppm) 


0.027 


0.216 


0.357 


0.373 


0.325 


0.287 


-0.649 


-1.422 


-0.876 


-1.848 



In the case of solution A, the relative frequency shifts Aff of A01 mode and A03 mode have similar 
properties with the solution of NaCl measured in Figure 3. But for the solution B and solution C, the 
relative frequency shifts of A03 mode is negative. As it is indicated in Equation (5), the absolute value 
pA$E of solution B and solution C measured in Table 2 should be smaller than the absolute value SsAp. 

Even for different kinds of solutions, such as the No. 6, No. 8, and No. 10 solutions, the values of Af/f 
of the A01 mode are almost the same, but the values of Af/f of the A 0 3 mode are apparently different. This 
means only one antisymmetric mode measurement, such as A01 mode, can not determine the solution 
with high selectivity. In order to specify the liquid with high selectivity, it is better to decouple the liquid 
density and the liquid sound velocity with the measured relative frequency shifts of A01 mode and A03 
mode. These will be analyzed in following section. 

3.4. Viscosity Measurements with the So Mode 

The responses of the So mode to the solutions of NaBr are similar to the ones of NaCl (Figure 3), 
the central frequency does not show apparent shifts due to their changing concentrations. However, the 
changes of the amplitude are related to the liquid's concentrations (Figure 4). It is because the 
frequency is mainly affected by the liquid density and sound velocity, which has no effects on the real 
part of the phase velocity of So mode (Equation (6)) [40]. The amplitude of the So mode in NaBr 
solution (Figure 4) shows that energy losses increase with the concentration, which is caused by the 
density and the viscosity of the liquid (Equation (7)). This can be used to decouple the density and the 
viscosity of the solution. 
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Figure 4. The responses of the So mode to different concentrations of NaBr solutions. 

-46.2 r , 




20.82 20.85 20.88 20.91 

Frequency (MHz) 

4. Results and Discussions 

4.1. To Decouple the Density and the Acoustic Sound Velocity with the Aoi Mode andAo3 Mode 

In order to specify the liquid with high selectivity, we attempted to decouple the liquid density and 
the liquid sound velocity with the measured frequency shifts of Aoi mode and A03 mode. Based on the 
Equation (4), two steps are undertaken in order to decouple the density and the sound velocity of 
a liquid as follows: 

(1) Determine the constant B { by measuring the frequency response of the reference liquid. 
Here, the water works as the reference liquid, the values of B\ and #3 are 14,767.3 N/m and 
1 0 1 ,642. 8 N/m respectively. 

(2) Get the physical parameters (pi, cl) for a unknown liquid by measuring cp,. All other 
parameters (B h A,-, M) in Equation (1) have already been calculated or measured. 

The first decoupled result using Equation (1) is the NaCl solution with different concentrations, as 
measured in Figure 3. Comparing the measured density and sound velocity with the values from the 
literature [32], the two results are close, as shown in Figure 5. 

Figure 5. The decoupled density and sound velocity of NaCl solutions and other solutions 
using measurements in Aoi mode and A03 mode. 
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Based on the same process, the density and the sound velocity for other three different species 
solutions with different concentrations were also decoupled and are shown in Figure 5. Compared with 
the values from the literature [32], No. 1-No. 6 solutions are KC1 solutions, No. 7 and No. 8 solutions are 
NaBr solutions, and No. 9 and No. 10 solutions are KBr solutions. Aside from decoupled density and 
sound velocity, the species of the solutions can be identified by comparing the measured values with the 
already known values [32]. After these density and sound velocity being decoupled using A 0 i mode and 
A03 mode, the liquid are determined with higher selectivity compared with only Aoi mode measurement. 

For Nos. 6, 8, and 10 solutions, the decoupled densities (Figure 5) are almost the same because the 
relative frequency shifts are very close (Table 2). For these solutions with adjacent density, sound 
velocities become the main factors affecting the frequency shifts in the A03 mode (Table 2). With the 
relations of the absolute frequency shifts: A/7/no.6 > A/7n 0 .8 > A/7/no.io (Table 2), the decoupled sound 
velocity of these three solutions have the following relation: Cno.6 > Cno.s > Cno.io (Figure 5). 

4.2. The Linear Response of (Density x Viscosity) 03 Changing with the Item of Amplitude Shifts (AA^b) 
in So Mode 

The item (density x viscosity) 0 ' 5 changes almost linearly with the item of amplitude shifts (A^b), 
and the linear fitting coefficient is about 1.02 + 0.09 (dB/(kg-m~ 2 -s~ 0 5 )), such as the NaBr and NaCl 
solutions (Figure 6). Beside the density and the viscosity, other parameters, such as the liquid's 
conductivity, affect the amplitude shifts. This is maybe the reason why the high linearity of 
(density x viscosity) 0 5 doesn't change with the amplitude shift. However, by analyzing the amplitude 
shifts (A^4dB) of the So mode, the viscosity of the solution can be derived. With the determined density 
by measuring the frequency response of the Aoi and A03 modes, the viscosity will be determined by 
checking the amplitude response in the So mode. 

Figure 6. Amplitude shifts versus (density x viscosity) 0 ' 5 for the NaCl and NaBr solutions 
in the So mode. 




(Densityx Viscosity) 0 5 (kg m" 2 s" 0 5 ) 

5. Conclusions 

In this study, for the first time to our knowledge, the density, sound velocity, and viscosity of the 
liquid are obtained simultaneously based on the measurements of the same solution sample volume 
with a Lamb wave sensor. When the frequency shifts of the Aoi and A03 modes are combined, the 
density and sound velocity are decoupled. This happens because the phase velocity of the A 0 i mode is 
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far from the sound velocity of the liquid and the phase velocity of the A 0 3 mode is close to the sound 
velocity of the liquid. Viscosity is obtained by measuring the amplitude of the So mode, which is 
especially suitable for Newtonian fluids. Effects of the aqueous electrolytic solutions' conductivity on 
these three modes are not clearly observed. The term (density x viscosity) 0 ' 5 doesn't show high 
linearity as a function of amplitude shifts. 

However, with this multi-parameter detection method, unknown solutions such as aqueous 
electrolytic solutions can be distinguished successfully and with high selectivity. Results clearly indicate 
that the multi-modes of a micro Lamb wave sensor are promising in the investigation of molecular 
thermodynamics, adiabatic compressibility, and molecular label free detection, among others. 
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